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Abstract 



X. 

H , Inclusive double differential multiplicities of deuterons, ^H, ^He, and ^He mea- 

sured by E877 for 11.5 A GeV/c Au+Au collisions at the AGS are presented. 
Light fragments at beam-rapidity are measured for the first time at AGS ener- 
gies. Beam rapidity deuteron and ^He yields and transverse slope parameters 
are found to be strongly dependent on impact parameter and the shape of 
the deuteron spectra is not consistent with that expected for a simple ther- 
mal distribution. The deuteron yields relative to proton yields are analyzed 
in terms of a simple coalescence model. While results indicate an increase 
in source size compared to collisions of lighter systems at the same energy, 
they are inconsistent with a simple coalescence model reflected by a rapidity 
dependence of the coalescence parameter B^. A new approach utilizing an ex- 
panding thermalized source combined with a coalescence code is developed for 
studying deuteron formation in heavy-ion collisions. The strong dependence 
of deuteron yields on collective motion implies that deuteron yields relative 



to those of protons can be used for constraining source parameters. 



Typeset using REVT^ 



INTRODUCTION 

The study of heavy-ion colhsions at ultra-relativistic energies has been pursued with the 
intent of observing matter at extreme temperatures and densities, the eventual goal being 
the production and study of deconfined and/or chirally restored matter. Although the 
relativistic heavy-ion program has developed over a number of years, and various hadronic 
observables have been studied, a complete characterization of the state formed during the 
high density phase of the collision is still outstanding. 

Measurements of light nuclei produced in the participant region of heavy-ion collisions 
are of interest since model calculations can be used to deduce measures of the source volume 
from the yield of nuclear fragments relative to that of protons |jl|-^. The yield of light nuclei 
is closely related to the two particle correlation between nucleons induced by final state 
interactions [Q . This can be incorporated into a theoretical framework from which one can 
deduce the density of the interaction region at freeze-out. 

In this publication we will report on investigations of the characteristics of spectra of 
light nuclei utilizing the E877 spectrometer, located at Brookhaven National Laboratory's 
(BNL) Alternating Gradient Synchroton (AGS). Spectra and yields of light nuclei, including 
deuterons, ^H, ^He, and ^He will be presented and compared to similar measurements at 
other beam energies and for various colliding systems. 

EXPERIMENTAL APPARATUS 

The E877 apparatus has already been described in [|]. Global observables, single parti- 
cle spectra of protons, pions and kaons [0 and two particle correlations HJ^ as well as 
collective flow [|10| have been studied by E877 at the AGS. The experimental setup primarily 



consists of three components: (1) a set of scintillators and silicon detectors for beam defini- 
tion; (2) a near Arc calorimetric coverage and a large acceptance charged particle multiplicity 
measurement, providing global event properties; and (3) a forward spectrometer, allowing 
measurements of identified charged particles. Figure [l| shows schematically the E877 exper- 
imental setup for the 1994 run. The z direction is defined to be along the beam and the y 
axis is defined to be out of the page (vertical) . 

The data presented here were collected in the fall of 1994. After event selection cuts 
were applied, approximately 5.4 ■ 10^ events for the 4% highest Et collisions were available 
for analysis. 

Beam Definition 

The beam in the C5 beam line at the AGS for the 1994 run had approximately 10'^ 
particles per spill with a Is spill length. To verify the beam composition and direction, and 
provide the start time of the experiment a number of detectors were placed upstream of the 
target. 

A series of four plastic scintillators (SI, S2, S3, S4) were used for the time reference 
definition in the experiment. Scintillators SI and S3 were ring-shaped veto detectors while 
scintillators S2 and S4 were eUipsoidal detectors placed at 45° relative to the beam. Fast rise 



time phototubes were utilized on S2 and S4 providing the interaction time with a typical 
Gaussian standard deviation of a = 35 — 40 ps. 

Two silicon microstrip detectors (BVER's) measured the beam angle in the xz-plane and 
were located between S2 and S4. The detectors were 300 /xm thick and composed of 320 
strips with a 50 /im pitch. This provided a measurement of the beam angle with a resolution 
of about 40 /xrad and determined the x coordinate of the beam at the target to 300 /xm. 
The BVERs were also used to tag and reject those events in which two or more particles 
traversed the experiment within the same 1 /xs time interval. 

To eliminate beam components with incorrect nuclear charge number, a silicon semicon- 
ductor counter (SILI) 87 microns thick was placed 2.5 cm upstream of the target to measure 
the energy loss of the beam particles. This detector provided a unit charge resolution of the 
incoming beam particles. 

Event Characterization 

For the 1994 run we used a AU target having a thickness corresponding to 1% of an 
interaction length. Two large, highly segmented calorimeters were used to quantify the 
centrality and reaction plane of the collisions. Together these detectors cover nearly 47r in 
the center-of-mass frame. 

The Target Calorimeter (TCAL) ||rT| consists of 832 Nal(Tl) crystals surrounding the 
target with a polar angle coverage of 48° < 6 < 135° (—0.88 < r] < 0.81). Each crystal has 
a depth of 13.8 cm, or about 5.4 radiation lengths (0.34 hadronic interaction lengths). 

The Participant Calorimeter (PCAL) |[ll| is a finely segmented lead/iron/scintillator 
calorimeter consisting of 16 azimuthal, 8 radial, and 4 depth sections. It covers a polar 
angle region of 1° < ^ < 47° (0.83 < 7] < 4.7). There is a small opening in the PCAL 
(—136 < ^2: < 16 mrad; —11 < 9y < 11 mrad) with iron wedges along the inside which 
define the acceptance of the forward spectrometer. 

Events are characterized by their centrality as determined by the transverse energy mea- 
sured by the calorimeters. 

Et = ^eisin6i (1) 

i 

where Cj is the energy deposited in one cell of the calorimeter and 6i is the polar angle of the 
cell with respect to the beam axis; the sum is performed over all PCAL cells. As was noted 
in [|11| , the overall systematic error associated with the measurement of Et was determined 
to be less than 4% for all collisions considered in this analysis. 

The transverse energy is then used as a measure of the centrality of the collision by 
noting that 

da s dN 



dEt UtB dEt 



(2) 



where a is the cross-section, B is the incident number of beam particles, s is the trigger 
down-scale factor [s > 1), and rit is the number of target atoms per unit cross-sectional 
area. For the experiment we used a Au target with an areal density of 540 mg/cm^. This 
corresponds to a 1% interaction probability for incident Au ions at an calculated inelastic 



cross section of (Tin = 6.1 b. Central collisions of a certain percentage (w) of the inelastic 
cross section were then selected by adjusting the lower limit E'""" such that [|l^: 

roo da ^771/ 
JSJ°™ dE'^-^t 

W = . (3) 

(Tin. 



Forward Spectrometer 

The momentum and time-of-flight of all charged particles accepted into the forward 
region are measured using two drift chambers, four multiwire proportional chambers and 
two time of flight hodoscope walls. These detectors are also shown schematically in Figure 

E- 

In addition, a set of multiwire proportional chambers was placed 2 m from the interaction 
region; the two chambers were separated from each other by 25cm (VTXA and VTXB). They 
were not used in this analysis. 

The analyzing magnet, situated downstream of the vertex chambers and spanning from 
z = 260 cm to 2; = 350 cm, could be run in two polarities in order to study possible 
systematic errors in the momentum measurement and to optimize the acceptance for positive 
and negative particles. The magnet had a maximum field of 3, 353 ±4 Gauss with a polarity 
that bent positive particles in the negative x direction and an integrated B ■ dL oi .3487 Tm 

i- 

Two drift chambers (DC2 and DC3), located at z = 700 cm and z = 1150 cm from 
the target, respectively, provided track x and y information and were used to determine the 
rigidity of each track; each drift chamber was composed of six planes of wires which measured 
the position in the bending plane of the magnet, i. e. the xz-plane, with a resolution of 250 /xm 
(350 /im) for DC2 (DC3). The anode wires in DC2 and DC3 were parallel to the ^direction, 
and the cathode planes in both chambers were segmented into chevron shaped pads. The 
signals from the pads provided y measurements of charged particles with a resolution of 
typically a few percent of the pad size. Each drift chamber had two segmentation size 
regions, resulting in ^/-resolutions of 2.3 mm (15 mm) in DC2 and 4.3 mm (36 mm) in DC3 
depending on the pad size. 

The multi-wire proportional chambers (MWPC's) with vertical anode wires space 5.08 
mm apart were situated between the two drift chambers. The MWPC's were used to ease 
tracking in the high multiplicity environment of Au+Au collisions. 

An array of 160 plastic scintillators (TOFU) |]13[ aligned vertically and placed behind 
DC3 was used to measure time-of-fiight with a time resolution in conjunction with the beam 
start counters of about 85 ps. It also provided an additional y-measurement with a precision 
of about 1 cm in the y-direction. The FSCI (Forward Scintillator) was a second hodoscope 
which also measured the time-of-fiight at a distance of 31 m from the target with a resolution 
of 350 ps. It was utilized in this analysis only to study of TOFU systematics since it had a 
lower resolution and smaller acceptance than the TOFU. 



DATA ANALYSIS 

Particle rigidity was determined from the radius of curvature of a track {R) through the 
magnet and the value of the magnetic field B: 

rigidity = ^^^^^ = RB (4) 

where Ze is the charge of the particle. The vertical momentum component, p^, is calculated 
from the vertical position of the track at the TOFU. The reconstructed mass of each particle 
is then given by 

where j3 and 7 are calculated from the time-of-flight measurement coupled with the calcu- 
lated path length to the associated TOFU scintillator. After all momentum calculations, 
Fig. 1^ shows the E877 acceptance for protons and deuterons as a function of transverse 
momentum and rapidity for the 1994 run. 

The charge of each particle is determined from the energy deposition in the TOFU slat. 
Shown in Figure ^ is the distribution of pulse-heights for those particles corresponding to 
the deuteron mass peak at beam rapidity. A clear Z = 1 peak is visible as well as a peak 
corresponding to Z = 2 ■ 1, from two Z = 1 tracks striking the same TOFU slat. In 
addition, the Z = 1 peak corresponding to ^He is clearly visible. The resulting distribution 
is fit to three Landau functions with a non-linear response due to the saturation of the 
photomultiplier tubes associated with the TOFU. The standard E877 analysis cut on the 
TOFU pulse-height distribution at 1.4 times the 1 MIP pulse height corresponding to about 
a 10% loss in the Z = \ yield with no contamination from Z = 2. 



Deuterons 

Since 5m/ m oc 5p/p and since the momentum and time of flight measurements worsens 
at larger momenta, deuteron selection is complicated by the tail of the proton peak at high 
momenta. To account for this effect, a systematic background subtraction technique was 
developed [|T^. In each rapidity and transverse momentum bin the yield of deuterons is 
calculated by fitting the mass distribution in this region to the proton and deuteron mass 
peaks plus a background between these peaks due to misidentified protons from resonance 
decays. 

Shown in Figure ^ is the reconstructed distribution of mass-squared in several pt = 
{px^ + Vy"^) and Udeut = 0.5\n{E + pz)/{E — pz) bins, where the rapidity was calculated 
assuming the deuteron mass. At low momenta, the deuteron and proton peaks are clearly 
distinguishable. In this region of phase space, the mass-squared distribution is described by 
two Gaussians and either an exponential or linear background^. At higher momenta, the 



^The systematic error associated with the choice of the background function has been found to 
be less than 5% . 



proton and deuteron mass-squared distributions broaden, due to the momentum resolution 
of the apparatus, and the background is minimal. In this region of phase-space, shown in 
Figure |B, the mass-squared distribution is described by two Gaussian distributions. In the 
2.6 < y < 3.0 region, the momentum resolution smearing coupled with the approximately 
equal particle yields for protons and deuterons, causes the peaks to overlap (Figure ^C), 
resulting in a loss of separation of deuterons from protons in this region. However, at 
deuteron rapidity y > 3.0 the contamination of protons is minimal because the proton 
would need to be traveling at twice the beam momentum. Thus, a deuteron peak is again 
identifiable (Figure ^D) and the yield can be determined with a single Gaussian. The widths 
of the proton and deuteron mass peaks were found to be consistent with the known time of 
flight and momentum resolution of the experiment. 

The yield of deuterons (A''^) in each small bin of phase space is determined from the 
parameters of a Gaussian fit to the deuteron peak: 



Nd = adAW27i/Am' (6) 

where a^^ and A^ are the variance and amplitude of the Gaussian distribution fit to the 
deuteron peak, and Am^ is the bin width. Those regions of momentum space that fall 
between the first two regions outlined above. Figure ^A and B, respectively, have been fit 
with both two Gaussians as well as two Gaussians plus a background in order to determine 
the systematic error due to the choice of fit funcion. The resulting yields of deuterons agree 
within 5%. 



Tritons 

The extraction of a triton signal is complicated by the fact that its mass-squared distri- 
bution is wider than the deuteron's due to multiple scattering and its yield is 1-2 orders of 
magnitude lower than that of the deuteron. However, a technique similar to that used for 
the deuterons can also be utilized to measure the yield of tritons in the collisions studied 
here. In Figure ^ is shown the reconstructed mass-squared distribution in various ytru bins 
for pt < 0.5 GeV/c. A clear triton peak is visible. The yield of tritons was calculated for 
each rapidity and pt bin by fitting this distribution with a Gaussian at the mass of the 
triton with an exponential background from the deuteron tail; the choice of an exponential 
or Gaussian background was found to have less than a 5% effect on the resulting yields. 
Several examples of these fits are also shown in Figure ^ 

Helium Isotopes 

Helium isotopes can be separated from Z = 1 species by cutting on energy deposition 
in the scintillators as described earlier. ^He are identified by recalculating the mass-squared 
according to eqn. ^ assuming Z = 2. This will cause all Z = 1 particles which pass the 
scintillator pulse-height cut to be found at double their actual mass. The ^He peak in this 
case is found between, and clearly differentiated from, the proton and deuteron residual 
peaks, respectively located at m"^ ~ 4 and m"^ ~ 14 (GeV/c)^, as shown in Figure ^. The 
^He peak is buried under the proton tail for y > 2.0 and, in the other rapidity bins, the 



yield of ^He may be extracted by fitting the mass-squared distribution by two Gaussians 
corresponding to the ^He and proton peaks. 

The method described above, i.e. cutting on the energy loss in the TOFU and, subse- 
quently, fitting the peak in the w? jZ'^ spectrum, was used to determine the deuteron, triton 
and ^He yield. However, this method cannot be used to determine the yield of ^He because 
the ^He peak in vn? jZ'^ space is buried under the significantly larger deuteron peak. We 
therefore use a complimentary approach to identify ^He: we first cut on the rn? /Z"^ peak 
and then fit the Z = 2 peak in the TOFU pulse height distribution. As described earlier, 
this distribution has two peaks, corresponding to one and two Z = 1 particles, at all but 
beam rapidity where an additional Z = 2 peak is evident as shown in Figure ^. These peaks 
are fit, at beam rapidity for pt < 2 GeV/c, to Landau distributions and the yield is then 
determined from the integral under the fit to the Z = 2 distribution. The overall systematic 
errors from this procedure is estimated to be less than 10%. 

RESULTS 
Mid-rapidity 

Shown in Fig. |^ are deuteron invariant multiplicities for the 4% highest Et collisions as a 
function oi pt and y with bin widths of 20 MeV in pt and 0.1 in y. The invariant multiplicity 
is fiat as a function of pt and y over most of the measured range. 

In Figure ^ are shown the measured invariant multiplicities of '^H and ^He for the same 
centrality. Plotted is the differential cross section at the center of the bin assuming an 
underlying thermal distribution with a shape similar to what is measured [|. The resulting 
measured cross section is not affected significantly by the assumed distribution. The distri- 
butions are presented as a function of rapidity for two pt regions. The error bars refiect the 
statistical uncertainty in the fit procedure and do not include the systematic errors, esti- 
mated to be about 20% for both ^He and ^H, dominated by the uncertainty in the assumed 
background. 

Beam-rapidity 

In Figure ^ are shown invariant multiplicities for deuterons at beam rapidity {ybeam = 3.2) 
for the 4% and 10-20% highest E^ events. As noted above, deuterons at beam rapidity are 
easily extracted from the mass-squared distribution. The transverse momentum distribution 
of beam rapidity deuterons for the 4% highest Et are distinctly harder, i.e. show a larger 
inverse slope parameter, than deuterons from the lower Et collisions. Over the range in pt 
in which we measure, a fit to an rrit Boltzmann distribution in the 3.0-3.1 rapidity bin gives 
Tb = 49.6 ± 1.4 MeV and Tg = 83 ± 2 MeV for the 4% and 10-20% Et bins, respectively. 



^The resulting slope and magnitude of the pt spectra was found to be insensitive to a wide range 
or reasonable assumptions of the true distribution in the acceptance correction. 



However, the reduced x^ of each of these fits is greater than 4, implying a Boltzniann shape 
does not describe the measured distributions well. 

In beam rapidity ^He spectra are shown in Figure |T0|. In the lower Et range (10-20%), 



the ^He yield is a factor of 10 higher than for the higher Et (4%), indicating that beam 
rapidity ^He are primarily produced by projectile fragmentation. 

MODEL COMPARISONS 

Light nuclei from heavy-ion collisions have historically been studied via a coalescence 
model [0,|15| which assumes cluster formation takes place at freeze-out. Since the binding 
energy of the deuteron is 2.24 MeV, it is easily broken apart in the fireball region where 



current models estimate the temperature to be 120-150 MeV [|I6],|T^. In the coalescence 



formulation [^ the probability of forming a deuteron is greatest when a proton and neutron 
at freeze-out have a small relative momentum. Assuming that the proton and neutron 
distributions are similar, the deuteron momentum distribution is then given by 

where the particle momenta obey pd = 2pp which can be easily generalized to heavier 
fragments with mass number A as 

^^^(P-) = B. [E.^iP.)) (8) 

where pa = App. In eqns. ^and ||, Ba {B^ in the case of the deuteron) is a phenomenological 
parameter known as the coalescence constant. In proton- nucleus collisions Ba has been 
related to the interplay between the binding energy of the deuteron and the optical potential 
of the target nucleus [jl5[. For nucleus- nucleus collisions it was recognized that the idea of 



a nuclear optical potential was no longer meaningful and Schwarzschild and Zupancic 
expressed the coalescence parameter in terms of a momentum difference pq. A proton and 
neutron would coalesce to form a deuteron if, at freeze-out, their momenta differ by less 
than pq. More generally, for the case of a source of Z protons and N neutrons with mass 
number of the cluster A, 



2^ / N\Z\ V 3 

where sa is the spin of the cluster and Rnp is the neutron to proton ratio in the source. 
Therefore, the coalescence parameter in these models is simply related to the momentum 
difference between corresponding protons and neutrons and is described by a simple step 
function in the coalescence probability at po- Note that no dependence on the colliding 
system exists in this model which described deuteron and ^H/^He production very well in 



nucleus- nucleus reactions at Bevalac |1T8| , P^ energies (about 0.1 — 2GeV-A) and in p- nucleus 
collisions, at FNAL [^ energies. For all of these energies, only a single constant Ba was 
needed to describe the light fragment distributions based on proton spectra in minimum bias 



collisions. However, for nucleus-nucleus collisions at AGS energies and above a dramatic 
drop below the Bevalac values in the coalescence parameter and a strong dependence on the 
centrality of the collision [§,0] was observed. For beam energies where nuclear collisions 
produce a large number of secondary particles, i. e. AGS energies and above, B^ was no longer 
independent of beam energy or composition presumably because the source decoupled at 
a size larger than the mean size of a deuteron. 



Rr 



2.1 fm f22|. In this case, B^ was 



theoretically deduced to be inversely proportional to the volume of the source at freeze-out. 
In the model of Sato and Yazaki [|l|, for example. 



B, 



3 (87r)3/2 

4 Z\N\ 



I'di' 



[i^d + i^) 



3/2 



(10) 



where Vd = -20 fm ^ is the Gaussian wave function parameter for a deuteron and u is directly 
related to the radius of the system at freeze-out by 



Rr 



Rx + Ry + Rz 



(3/2Z.: 



,1/2 



(11) 



where Rx, Ry and R^ are the Gaussian radii. Other models [Q also utilize a wave-function 
description and a density matrix formalism to describe deuteron formation. With the ex- 



ception of large cascade codes p2|j2^ as well as recent hydrodynamically motivated source 
parametrizations [^, however, none of them explictely account for transverse collective 
expansion. Though such an approach is not excluded from the Sato and Yazaki approach 
|l|], nothing beyond sources with no space-momentum correlations were considered in their 
presentation. Further, though there have been some attempts at a fully relativistic treat- 

this generalization is non-trivial because of difficulties in the proper 



ment of coalescence 26 



relativistic treatment of bound states 



Coalescence Parameter Ba 



The coalescence parameter, as defined in Eq. 0, can be calculated from E877 measure- 
ments of protons P and deuterons. To properly calculate B^ for our data set, the fraction 
of protons from A decays must be subtracted from the measured proton spectra in order 
to obtain the proton distribution at freeze-out. This correction was performed utilizing the 



measured A yield in Au-|-Au collisions |2^ and knowledge of the acceptance and reconstruc- 
tion efficiency of the experimental apparatus. Over the measured acceptance region, the 
hyperon contribution to the proton spectra was found to be no more than 10%. 

After correcting for hyperon contributions, Bd can be calculated as a function of pt and 
y. We notice, within the acceptance of the present measurement, no variation of Bd as a 
function oi Pf However, there is a significant rapidity dependence of the average Bd for all 



measured pt as shown in Figure |rT[ If Bd is related to volume at freeze-out, such a result 
may imply a changing effective freeze-out volume as a function of rapidity. The measured 
change in Bd of nearly a factor of two from mid-rapidity to y = 2.6 would correspond to a 
decrease in the effective source radius of 1.25. 

Figure |T^ shows a comparison of the coalescence parameter deduced from experiments at 
Bevalac [||,||] energies with those obtained at the AGS in in Si + Al jSH, Si + Pb p8|, Au + 



10 



Pt p9[, and Au + Au |3^ collisions. The error on the Au+Au data point corresponds to the 



minimum and maximum range of B^ values shown in Figure |ll] including the statistical error 
bars. All presently measured values of 5^ for high Et Au+Au collisions are significantly lower 
than the results for Si+Al and Si+Pb. This is consistent with the interpretation that B^ 
depends on the volume at freeze-out. Two proton and two pion correlations measurements 
have found that the effective nucleon source size is larger for the Au+Au system than for 
lighter systems 0,^]. 

The coalescence parameter measured in these collisions is nearly one-tenth of that mea- 
sured in central Au+Au collisions at Bevalac, implying in the models of P],p3[ that the 
volume at feeze-out in AGS Au+Au collisions is five to ten times that at Bevalac energies. 
Such a simple geometrical interpretation would be inconsistent with two particle correlations 
measured at the AGS ^^■ 

A possible reason for such an inconsistency is that the models which deduce an inverse 
volume dependence [0)0 ^'^ ^°^ account for collective motions in the source. As has been 



shown in [^, the yield of deuterons is highly dependent on the amount of collective motion 



in the source. Hence, B^ is not simply related to the volume but also the position-momentum 
correlations at freeze-out. If the source is perfectly correlated, then proton-neutron pairs 
that freeze-out near one another in space will have similar momentum vectors. The deuteron 
yield would be larger in this case than for a random source. As such, any consistent model 
of deuteron production must account for correlations in the source. 

Therefore, in order to completely describe deuteron production, a more sophisticated 
model of the source that includes dynamical correlations must be formulated. In large 
cascade codes it is difficult to study the influence of the source parameters. We thus used a 
simple dynamical model in the following section to study such effects. 

In Fig. O is the corresponding coalescence parameter for ^H and ^He extracted from this 



data set compared with Bevalac [|T8[, AGS [g,0, and SPS |31|] values. Similar to what is 
observed in Bd, the three nucleon coalescence parameter drops dramatically at AGS energies, 
consistent with an increase in the effective source size. 



A DYNAMICAL MODEL FOR DEUTERON FORMATION 

An analytic, hydrodynamically motivated model, devised by Chapman, et al. pH] has 



been implemented to describe the particle source distribution in heavy ion collisions [34 



This model uses an emission function (5*) for a flnite, expanding, locally thermalized source 

Ek^^ = jd'xS{x,K) (12) 

where x and K are the position and momentum 4-vectors, respectively. The source includes 
longitudinal and transverse collective motion with an underlying temperature, given by 



Six K) - ^^^"^M^-^) ,^, 



K ■ u{x) (r - ToY r^ {t] - t]o) 



2(Ar)2 2i?2 2{Ar]y 



:i3) 



where 



11 



V = \H{t + z)/{t-z)] 



T 



^{t' - z') 



r = ^/n[(i + A)/(i-A)] 

X 

K ■ uix) = Mt coshfr] — Y) coshrjAr) — Kt— sinh r^^fr) 

r 

Vt{r) = rif{r/R). 

In this model (x, y, z, t) is the freeze-out position four vector, r = ^/x^^^P^p and /?; and 
Mf are the longitudinal velocity and transverse mass, respectively. Kt is the transverse 
momentum and rjf is the transverse flow velocity. Furthermore, Ar, A?7 and R describe the 
respective Gaussian widths of the proper time, longitudinal velocity and one dimensional 
Gaussian radius. For more information on the source function, see p3| , |3^ . 

To determine deuteron yields with this framework, a program devised by R. Mattiello 



3^ , |36[] was implemented. It utilizes a method that has been applied to bombarding energies 

from 1 GeV per nucleon in association with the the intranuclear cascade model and QMD 

3^ up to AGS and SPS energies with the cascade models RQMD [§^,0 and ARC g. In 



this approach, the number of deuterons is given by a summation over all proton and neutron 
pairs at freeze-out accounting for the Wigner density of the Hulthen deuteron wave function 



2^ . This method is equivalent to the density matrix formalism of Sato and Yazaki [|1| in 
the case of a static source. Therefore, by projecting the neutron/proton position-momentum 
distribution on the deuteron wave-function via this Wigner method, we may determine the 
distribution of deuterons in phase-space given any particular source of nucleons. 

Shown in Figure 113 is the distribution of the reduced x^ minus the minimum x^ {x^m = 



0.96) with respect to the measured deuteron yields calculated as a function of T and R with 
a resolution of 10 MeV and 0.25 fm respectively. For this study we assume that tq = 3fm/c 
and Ar = Ifm/c in Eqn. |13|. [| The technique we utilize for mapping out this x^ space 
is the following: For each value of temperature T and one dimensional Gaussian radius 
i?, the transverse and longitudinal proton spectra uniquely define the transverse (t^j) and 
longitudinal (Ar^) boost velocities. Every value of T and R in Figure |14| can describe the 
proton transverse momentum spectra given a particular rjj at mid-rapidity. 

We have found that the proton rapidity distribution in this model poorly reproduces 
the data [|5|,P5|, with a reduced x^ oii the order of 10 for all choices of T, R and Ar^. 
Furthermore, the shape of the transverse momentum distribution of nucleons in this model 
is assumed to be constant as a function of rapidity, in conflict with experimental data [Q . 
The assumptions inherent in the model I^S] regarding a hydrodynamical source are not 
entirely valid for nucleons at AGS energies due to these complications in the longitudinal 
direction. Therefore, in the studies presented here we only consider the nucleon production 
at mid rapidity and implement this model by fitting rjf to the proton pt distribution only at 
mid-rapidity and constrain At] to approximate the rapidity distribution. The deficiencies in 



^For SPS energies, Ar ~ 1.5 — 3.0 fm/c from two particle correlations [p^j35| while ||2^ noted 
that tq ~ 9fm/c. Similar numbers for AGS energies have not been published. 
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the model near beam-rapidity should not influence the results we show here for mid-rapidity 
and indeed the model describes the transverse momentum distribution at mid-rapidity quite 
well (x^ ?a 1). 



For this analysis, the proton spectra from E866 Au+Au collisions ||3g] are utilized due 
to their coverage of protons at mid-rapidity. The model is therefore fit with these two 
parameters to the proton spectra. After the proton distribution is determined, the deuteron 
yields at mid-rapidity are calculated with the coalescence algorithm described above. Results 
are compared to deuterons measured by E877 at mid-rapidity and pt < 0.6 GeV/c. 

The resulting x^ distribution strongly constrains the temperature and transverse radius 
in this model of the nucleon source though the two variables are highly correlated. Studies 
of proton correlations for the Au+Au system at AGS energies p| concluded there was an 
increase of the transverse freeze-out distribution of nucleons beyond the Au rms radius of 
S.lfm. This is consistent with the current results where all reasonable values of temperature 
result in a value of the transverse radius significantly larger than the rms radius of a Au 
nucleus. Furthermore, measures of the temperature at freeze-out in Au+Au collisions at the 



AGS |T^,|T^ of 120-150 MeV are also consistent with the results of this study. 

Note, however, that this method has not constrained Ar or tq. Varying these attributes 
from the values chosen above by a factor of two can vary the yield of deuterons by 10%. 
An increase in either Ar or tq would decrease the yield of deuterons in the model and, as 
a result, would shift the x^ curve of Figure O to lower temperatures and radii. Therefore, 



deuterons cannot be used to exclusively determine source parameters but must be used in 
conjunction with complimentary observables to produce a coherent view of the source. 

CONCLUSIONS 

We have measured light fragment yields at both mid and beam rapidities. Invariant 
multiphcities have been presented for deuterons, ^H, and ^He at mid-rapidity. Deuterons at 
mid-rapidity have been interpreted in a coalescence framework. The resulting coalescence 
parameter (5^) is consistent with an increase in source size from collisions of lighter nuclei 
at the same energy. However, a clear dependence on rapidity of -B^ implies the assumption 
of a constant effective volume as a function of rapidity is incorrect. 

We have introduced a technique for using deuterons as a sensitive constraint on source 
parameters in a heavy ion collision by incorporating a simple, hydrodynamically motivated 
nucleon source function at mid-rapidity coupled with a coalescence algorithm. The results 
show that deuteron yields are sensitive to a convolution of volume and collective effects, 
similar to other correlations measurements. These results can be used in conjunction with 
other measures of the source to tightly constrain theoretical models and assumptions. 

First measurements of beam rapidity deuterons and ^He at AGS energies have also been 
made. The resulting distributions have a strong Et dependence, consistent with a production 
mechanism dominated by spectator fragmentation. 
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the 1994 run. 



CZ3 

•I— I 

i 

•1— I 




2 2.5 3 3.5 4 

Normalized Pulse Height 
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